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Abstract. 

l/^ , We present spectra of twelve Type la supernovae obtained in 1999 at the William Herschel Telescope and the Nordic Optical 

Telescope during a search for Type la supernovae (SN la) at intermediate redshift. The spectra range from z = 0.178 to 
z = 0.493, including five high signal-to-noise ratio SN la spectra in the still largely unexplored range 0. 15 < z < 0.3. Most of 
the spectra were obtained before or around restframe B-band maximum light. None of them shows the peculiar spectral features 
found in low-redshift over- or under-luminous SN la. Expansion velocities of characteristic spectral absorption features such as 
Sill at 6355 A, SII at 5640 A and Call at 3945 A are found consistent with their low-z SN la counterparts. 

£: 
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1. Introduction Although the number of high-redshift SN la discovered 
' and followed-up greatly increased along the years and is now 
d ■ Type la supernovae (SN la) provide us with a powerful tool rapidly approaching a ^^1, the gap observed in the 
for constraining cosmology through the magnitude-redshift Hubble diagram at intermediate redshift range has endurecL 
Hubble diagram. In the last decade, systematic searches for i ndee d, optimizing supernova searches in that redshift range re- 
SN la have been performed in view of increasing the statistics quiref , securing , arge fractions of observing time on an interme . 
at all redshifts, while at the same time trying to improve our diat£ size telescope ^ 2m) with a large field of view ( > 2502) 
understanding of these events. Programs launched in the early in order tQ ensure effidmt detection of a large enough num . 
nineties were targeting high-redshift (0.4 < z < 0.8) supernovae ber of supernovae for the duration of the survey . intermediate 
and have lead to the discovery of the acceleration of the expan- redshift supernova programs are thus subject to the availabil- 
sion of the universe j Riess et alJ | 1998HPerlmutter et alj | 1999|) . ky of such telescopes and instrume nts. Nevertheless, observ- 
More recently, even more distant SN la (up to redshift ~ 1.7) ing supernovae at interm ediate redshift is crucial not only to 
have been discovered and followed-up with HST/ACS allow- „ bridge ±Q gap „ between ^ w _ and high . redshift samples 
ing one to "observe" the turn-around epoch of the acceleration tQ increase ±e precision in measuring the cosmological pa- 
history of the universe j Riess et aU | 2004|) . The nearby super- rameterSi but also t0 obtain precise rest f rame U-band photom- 
novae (z < 0.1) have also received great attention, as they are gtry and spectroscopy of SN Ia> otherw ise difficult to obtain on 
necessary to anchor the Hubble diagram at low redshift and be- 

cause they offer the possibility to study systematics since they ' current high-statistics programs such as ESSENCE 
are easier to observe than their high redshift counterparts. (www.ctio.noao.edu/~wsen) or SNLS (cfht.hawaii.edu/SNLS) 
now target number statistics of well measured SN la well above the 



Send offprint requests to: christophe.balland@ias.u-psud.fr several hundreds in large multi-year projects 



2 



Balland et al.: Spectroscopy of twelve Type la supernovae at intermediate redshifts 



nearby objects because of the low detector efficiency and strong 
variability of the atmosphere in the U-band. 

Of utmost importance is to check that the properties of 
SN la do not significantly evolve with redshift. Among po- 
tential systematics is the evolution of SNe la spectral char- 
acteristics with z. Differences in the SNe la environment and 
parent population are likely to translate into differences in 
the presence and evolution of spectral features. Indeed, as 
our knowledge of observable properties of SN la improves, it 
clearly appears that diversity rather th an uniformity is common 
among type la supernovae (|Branchl|l987t [B ranch et all Il988l 
lHatano et alJl2OO0t iBenetti et alJl2004 120051) . It is crucial to 
quantify this diversity in order to estimate to what extent it af- 
fects the validity and accuracy of the cosmic parameter deter- 
mination from the Hubble diagram. In this respect, an enlarged 
sample of precisely measured intermediate redshift SN la spec- 
tra is needed to improve our understanding of their behavior. 

It is thus important to discover and follow-up both pho- 
tometrically and with spectroscopy a substantial number of 
SN la in the intermediate (0.1 < z < 0.4) redshift range. 
For this reason we launched, during the year of 1999, a cam- 
paign to search and follow-up SN la in that redshift range. A 
second campaign was pursued during the year 2002. In this 
paper, we present the spectra taken at the William Herschel 
Telescope (WHT) and at the Nordic Optical telescope (NOT) of 
the twelve supernovae discovered in the spring and fall 1999, 
at the Isaac Newton Telescope in La Palma observatory. The 
photometric follow-up of five of these supernovae and their in- 
clusion in the Hubble diagram will be presented in a subse- 
quent paper. In Sect. 2, we briefly describe the search strategy 
and discovery technique used. In Sect. 3, the spectroscopic ob- 
servations of the supernova candidates, including the reduction 
procedure, are presented. In Sect. 4 and 5, we focus on the iden- 
tification and analysis of our set of SN la spectra. Discussion 
of results and comparison to other works are done in Sect. 6. 
Concluding remarks are given in Sect. 7. 

2. Search for SN la at intermediate redshift 

As the present paper focuses on the spectral analysis of the 
SN la sample, we only briefly describe here the strategy and 
technique used to discover the supernovae. 

2.1. Overall strategy 

The search for supernovae was performed on images taken us- 
ing the Wide Field Camera (WFC), a four 2kx4k thinned EEV 
CCD mosaic with a ~ 0.27 square degree FOV mounted on the 
Isaac Newton telescope. Most of the observations were done as 
part of the Isaac Newton Group's Wide Field Survey (WFS), 
a public access survey of multicolor data obtained over a 5 
year period through the INT photometric filters set (u' through 
z') over more than 200 square degrees with a typical d epth 
of r' * 24 and g' * 25 JWalton etal.ll 19991 iMcMaho netail 
2001 ). WFS observations were scheduled as to permit repeated 
observations of the same area of sky. Detection of the SN candi- 
dates was done by repeated imaging, in the Sloan-Gunn g' 1996 
filter, of the South Galactic Cap (SGC) fields over an area of 



-31 square degrees. SGC fields at zero declination were cho- 
sen so that follow-up from both northern and southern hemi- 
spheres could be done. After a test period in April 1999 whose 
purpose was to validate our strategy and detection technique, 
600 s images were taken around the new moon in August 1999 
(new moon on August 1 1th) and then used as the input baseline 
reference images for the subsequent supernova search. As the 
restframe rise time for S N la is about 18 days dAldering et alJ 
2000; Ri ess et al.ll999al) . two new images of 240 s each of the 
same fields in the same g' filter were taken roughly a month 
later (around the September new moon) in order to discover 
the s upernovae around their maximum light ( Perlm utter et alJ 
119971) . The schedule of this new set of observations was opti- 
mized to ensure the best matching in seeing between the refer- 
ence and new images. The image quality of the data was gen- 
erally good, with a median seeing of as 1 arcsec. 



2.2. Discovery technique 



The search procedure generates potential SN la candidates, as 
well as other objects, such as AGN, other SNe types, asteroids, 
etc. The new image (the discovery image) consisted in two ex- 
posures taken one hour apart, so as to reject cosmic ray hits and 
asteroids. After pre-processing , the discovery images were co- 
added and three subtracted frames were built, by subtracting, 
from the reference image, one or the other discovery image, 
or their sum. These subtracted frames were then searched for 
point-like sources. The supernova candidates were selected by 
requiring that they were detected on the 3 subtracted frames, at 
the same place within 2 pixels (0.66 arcsec) in order to reject 
cosmic rays and asteroids, and that their flux within an aper- 
ture of 1 full width at half maximum (FWHM) in radius is at 
least larger than 15 % of the flux measured in a similar aperture 
on the reference frame (so as to eliminate subtraction resid- 
uals). This requirement also rejects slow varying objects, as 
plateau type II supernovae. Finally, the supernova candidates 
were visually checked by two independent scanners. These 
candidates were re-observed the following night, and kept only 
if re-detected with a similar (or higher) flux. This last proce- 
dure tends to eliminate any slow-moving object and supernovae 
on the decline. The quality of the candidate as a potential su- 
pernova was estimated from its flux and shape parameters and 
prioritized as very likely, probable or possible supernova. A 
total o f fifteen candidates wa s obtained (the 'September candi- 
dates' lAstier & Goobarl Jl999l) '). Spectra were obtained within 
a few days after detection during pre-scheduled spectroscopic 
follow-up time at WHT. For two candidates, complementary 
spectra were obtained at the Nordic Optical Telescope (NOT) 
about three weeks after discovery. Ten out of the fifteen spec- 
troed candidates have been confirmed as SN la . To this set, we 
add two more supern ovae (the 'April candidates' IHardin et alJ 
lll999h :l Hardin ( 1999)) obtained in a similar way during the test 
run in April 1999 (see below). 
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3. Spectral observations and reduction 

3.1. Observations 

Spectra of all the supernova candidates were obtained with 
the ISIS spectrograph at the 4.2-m William Herschel telescope 
(WHT) on April 19 and 20 (SN 1999cj and SN 1999ck discov- 
ered during the test run) and on September 12th, 13th and 14th 
1999 for the fifteen September candidates. 

For April 1999 candidates, both the red and blue channels 
of ISIS were used along with the ISIS 6100 dichroic filter. For 
the blue part, we had a blue sensitive thinned EEV CCD of 
4096x2048 pixels combined with the grating R158B (useful 
spectral range from 4000 to 6200 A, dispersion 1.62 A per 
pixel). For the red part, the detector used was a 1024 x 1024 
pixel TEK CCD device operating in its standard mode and 
combined with the R158R grating covering a total useful spec- 
tral range of ~ 6000 to 8000A with a wavelength bin per pixel 
of 2.90 A. With this device, fringing is negligible below 7500 
A. 

Three exposures of 1800s each were obtained for the two 
April candidates, both in the blue and red parts, with a long 
slit of 1 arcsec width. Copper-neon arcs were used for wave- 
length calibration. Two standards were observed and used for 
flux calibration (BD+33°2642 and BD+28°421 1, see below). 

During the September campaign, only the red channel of 
the spectrograph was used in order to avoid some sensitivity 
loss due to the dichroic filter. Even though the two April spectra 
were of acceptable quality, removing the dichroic filter allowed 
us to optimize the efficiency in the red channel, where the prin- 
cipal features of SN la are present for our redshift range. The 
red arm had the same detector as in April. Three long-slit spec- 
tra of 1800 s each were obtained for each candidate, except for 
SN 1999dx, SN 1999dv and SN 1999gx for which only two 
1800 s exposures were taken. A slit width of ~1.5 arcsec was 
used. Parallactic angle was generally used to minimize loss by 
differential refraction. In some cases, the slit was put at a spe- 
cific orientation in an attempt to take a spectrum of the host 
galaxy at the same time as the supernova spectrum. Three spec- 
trophotometric standards were observed each night (Feige 25, 
BD+28°421 1 and Feige 1 10), with both -1.5 arcsec and 8 arc- 
sec widths. Copper-neon arcs were obtained at the beginning 
and end of each night for the purpose of wavelength calibra- 
tion. 

Weather conditions were good for the three September 
nights with a seeing of about 1 arcsec or less, except for the first 
half of the second night (seeing about 3 arcsec). Five candidates 
were observed during each night. Two additional exposures of 
one faint candidate (SN 1999dz first observed on September 
13) were obtained during the night of September 14th. 

Additional observations of SN 1999du and of SN 1999dv 
were obtained on October 2nd, 1999 at the 2.5m Nordic Optical 
Telescope (NOT) in La Palma. Long slit spectra have been ac- 
quired with the ALFOSC spectrometer, using a 300 lines/mm 
grism (number 4) and a CCD Loral (2048x2048 pixels) detec- 
tor of 15 fxm pixel size. The wavelength coverage ranges from 
4200 to 8900 A. However, due to the presence of the second 
order spectrum with this grism, the useful range is reduced to 



approximately 4200 to 6000A. A slit width of 2.5 arcsec and 
of 1 arcsec, yielding a resolution of 35-40 A and of 14 A have 
been used for SN 1999du and SN 1999dv respectively. Five 
spectra of SN 1999du and six spectra of SN 1999dv were ob- 
tained with an exposure time of 1800 s. The last spectrum of 
SN 1999du is very weak and not considered in the following. 
BD+28°4211, observed on October 3, has been used for flux 
calibration purposes and helium lamp spectra were taken for 
wavelength calibration. The seeing during the NOT observa- 
tions was always better than 1 arcsec. 

A summary of the spectroscopic observations are given in 
Table[T] 

3.2. Reduction 

A preliminary assessment of the spectroscopic data was per- 
formed interactively at the telescope to allow for confirmation 
of the candidate, so as to select the supernova sample for pho- 
tometric follow-up as soon as possible after discovery. 

The spectroscopic reduction was redone using the ESO- 
MIDAS data reduction software version 02SEP. An average 
bias and normalized flat field frames, obtained from an internal 
tungsten lamp, are produced for each night using a median filter 
in order to remove cosmic ray hits. A two dimensional wave- 
length calibration is done, using copper-neon arcs, by comput- 
ing a dispersion relation for each row of pixels along the disper- 
sion direction (horizontal with our set-up) over the entire two 
dimensional frame. This corrects for some distortion along the 
wavelength direction and minimizes subtraction residuals after 
sky removal. About 40 lines are identified and used, covering 
the total spectral range in a regular way. This leads to a pre- 
cision better than 0.05% over the whole spectrum. Finally, we 
check on a few prominent sky emission lines (Nal at 5893A 
and [OI] at 5577, 6300 and 6363 A) that the wavelength accu- 
racy is consistent with this value, corresponding to less than 1 
A difference between the calibrated sky line and its theoret- 
ical wavelength. Extraction of the spectrum is then performed 
from th e 2D calibrated frame using a procedure based on lHornel 
( 1986) weighting algorithm optimizing sky subtraction and flux 
extraction. An error spectrum, essentially dominated by the sky 
signal, is generated. 

After correcting for atmospheric extin ction using the ex- 
tinction curve of La Palma observatory dKindfl985l) . abso- 
lute flux calibration is done. A response curve is computed 
by comparing our observed spectrophotometric standard stars 
with their tabulated absolute flux (in practice, we used Feige 
110 to calibrate all the September spectra, and BD+28°4211 
and BD+33°2642 for, respectively, the blue and red part of the 
April spectra). This response curve is then applied to our spec- 
tra. Finally, we remove "by hand" sky residuals resulting from 
bad sky subtraction or remaining cosmics. We do not attempt 
to remove atmospheric absorption lines around 6900 A. 

For the two April spectra, the reduction and extraction was 
performed on the blue and red frames in a separate way. After 
flux calibration, the two parts of the spectrum were combined to 
give the final full spectrum, with a useful wavelength coverage 
from about 4000 to 7500 A. The two parts overlap nicely over 
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~ 200 A around 6100 A , showing that our flux calibration or 
extraction procedure is not significantly in error. 

As for the WHT data, standard reduction for NOT data was 
performed in the ESO-MIDAS environment (cosmic removal, 
bias subtraction, flat-field correction, 2D-wavelength calibra- 
tion, and flux calibration). Again, the wavelength calibration 
was checked against sky lines. They were found within 3-4 A 
from their tabulated wavelengths. Extraction was done follow- 
ing the same procedure as for WHT spectra. 

4. Identification 

4.1. Redshift determination 

Redshift determination is based on the identification of charac- 
teristic galactic emission and absorption lines which are nar- 
rower than any supernova feature. We use mainly [Oil] at 
3727A, Call H&K at 3945 A, [OIII] doublet at 5000A or Hp, 
or several of these lines together if present. We do a gaussian 
fit of the line profile. The relative spectral resolution 6A b s /A b s 
is typically less than 0.05 % over the useful spectral range (see 
Sect. 3.2) which translates into an absolute error on the redshift 
6z - SA bs/ AobsXil +z)less than 0.001 atz = 0.5 and less than 
0.0005 at z ~ 0.1 . In the following, redshifts are thus given with 
three significant figures. 

4.2. SN template fitting 

Identification of the candida tes is done using a software tool 
developed by one of us ( Sainton 2004a b). The software, called 
SN-fit, is based on a x 2 minimization fitting of the reduced 
spectrum to a set of known spectral templates. As contamina- 
tion from the host galaxy often occurs, a model spectrum M 
consisting of the weighted sum of a fraction of a supernova S 
and a fraction of a galaxy template Q (or of the observed galaxy 
host when available) is built in the restframe: 

M(A resr ,z,a,/3) = aS(A rest (l+z)) + ffi^restV + z)) (D 

All template spectra are initially in the restframe and are appro- 
priately redshifted. The redshift determined from galaxy lines 
(see previous subsection) is used to constrain the input redshift 
interval to search for the solution, although it is not manda- 
tory as the input redshift interval can be left arbitrarily large (z 
is then considered as a free parameter). The x 2 fitting is done 
against the set of parameters (z,a,/3) and is robustified against 
undesirable points such as regions where atmospheric absorp- 
tion takes place or where sky subtraction residuals are present. 

The fit is done for every template couple (supernova, 
galaxy) available in our database and results are classified in in- 
creasing order of reduced^- 2 . In order to assess the significance 
of the best-fit solution, we systematically examine the imme- 
diately following solutions to check for possible degeneracies 
in the models. It is clear that the result will greatly depend on 
the content of the database and special care has been put to 
build it. It contains about 150 spectra (mostly supernovae plus 
a dozen galaxy templates) taken from various sources. The set 
of supernova templates contains local 'normal' SN la as well 
as several peculiar la and Ib/Ic/II supernovae. These spectra 



have been put into the r estframe. Galaxy templates are taken 
from Kin nev et alJ ( 1996). Morphological types include bulge, 
E, SO, Sa, Sb, Sc, and starbursts (Stb) with various amount of 
reddening JCalzetti etalJll994 iKinnev et ai1ll996h . We refer 
the reader to Appendix A for more details. 

One concern with the database is the uniformity of the sam- 
pling of the supernova phase (the date of the SN spectrum 
template with respect to maximum light). If the sampling is 
too scarce, this limits the resolution on the phase of the ob- 
served spectrum we get from the fit (we call this phase the 
"spectral phase"). We find a typical +3 day difference on the 
spectral phase. Higher differences between two solutions sim- 
ilar in terms of x 2 can happen when the SN signal is weak 
with respect to the host galaxy or when the spectrum of the 
SN is acquired well after maximum light, as spectral features 
evolve less rapidly a few weeks after maximum (phase degen- 
erate solutions). Since the detection procedure is optimized to 
pick up SN around maximum, this latter case happens only 
rarely. Another concern is related to the fact that we do not try 
to make any accurate determination of the reddening to the SN. 
Indeed many effects combine to account for the reddening that 
are difficult to separate: differential slit losses when not at the 
parallactic angle, errors in flux calibration and template spectra 
uncorrected for reddening. Mode ling 'reddeni ng' by adding a 
term in equation (1), as in lHowelletalJfc002ll .' is thus subject 
to caution as interpretation of such a term is difficult and we 
choose not to do so. This means that some degeneracy might 
exist between the best-fit phase and a possibly 'reddened' SN 
or template. This might increase the uncertainty in the phase 
determination. 

To quantify this effect, we have artifici ally reddened th e 
template spectrum of SN 1994D at - 2 days JPatat et aljH~9 96). 
using the reddening law of Howarth jHowarthl ll983) for E(B- 
V)=0.2, 0.4, 0.6 and 0.8. We have then fitted this spec- 
trum, usi ng .SAT-fit, w i th a s et of spectral tem plates intro- 
duced by iNugent et al.l ll2002h and modified by iNobili et alJ 
(2003). These templates are constructed from a broad sample 
of SN la and they are calibrated so as to reproduce the mag- 
nitude evolution along the supernova light-curve in the BVRI 
passbands. Spectrally, they offer, contrary to local templates 
of individual SN la , a phase resolution of 1 day over a large 
phase range. For our purpose, we only consider templates be- 
tween -10 and + 10 days. For each value of E(B-V), we com- 
pute the phase from the first best-fit solutions by weighting the 
phase 4>i of each solution with the corresponding x 1 probabil- 
ity For the unreddened spectrum of SN 1994D, as well as 
for E(B-V)=0.2, the phase obtained is « -3.5 days. The orig- 
inal -2 days phase is thus recovered within less than 2 days. 
For E(B-V)=0.4, 0.6 and 0.8, the difference between the origi- 
nal phase and the determination from SN-fii increases: the ob- 
tained phase is m +3.5, +5.5 and +9 respectively. As the origi- 
nal spectrum is progressively reddened, later-phase (i.e. redder) 
spectra are selected by the identification procedure. Note that 
the visual aspect of the fit is considerably degraded for E(B-V) 
values higher than 0.2. We conclude that for moderate redden- 
ing, the phase determination is reliable within 2 to 3 days. For 
higher values, no satisfying fit is obtained. In these cases, vi- 
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sual examination of the spectrum reveals the reddened nature 
of the spectrum. 

An other concern is the possible misidentification of a 
SN la spectrum around maximum as a type Ib/Ic or II super- 
nova. The P-Cygni features of a type II spectrum might fit some 
portion of the red part of a SN la spectrum before or close to 
maximum. More likely, a type Ic spectrum might be misiden- 
tified as a type la. In particular, this can happen when the fit is 
performed on a too small portion of the spectrum to be identi- 
fied. In any case, it is important to have a sufficient number of 
type Ib/Ic/II spectral templates in the database. At the time of 
this analysis, we dispose of 64 spectra of type Ib/Ic/II objects 
(see Appendix A). 

In the present analysis, the redshift is always determined 
from galaxy lines and it can be safely constrained in SN-fit. 
Templates selected in the spectral database typically cover a 
large wavelength range but as they are redshifted, the overlap 
between the model spectrum (redshifted SN + Galaxy) and the 
observed spectrum might be small. This can lead to an unreli- 
able identification as confusion with other types may occur. It 
is thus important to check a posteriori that the best-fit solution 
spans a large enough wavelength range. In case it is too small, 
inspection of adjacent solutions is necessary to assess the va- 
lidity of the best-fit solution. 

Obviously, the reliability of the identification depends on 
the underlying model, which is directly connected to the num- 
ber, quality and variety of templates available in the database. 
A good phase sampling is necessary but a sufficient number of 
spectra of different SN la supernovae at a given phase is also 
required in order to reflect the intrinsic differences observed in 
SN la spectra at low redshift. Note that the higher the quality 
of the observed spectrum is, the more significant the discrep- 
ancies between the model and the observed spectrum are, in 
terms of x 2 value. A high best-fit x 2 value for a given spec- 
trum either happens when the observed spectrum does not cor- 
respond to a supernova (and thus no satisfying model can be 
constructed from templates available in the database), or when 
the signal-to-noise of the observed spectrum is high. In any 
case, x 1 changes between solutions for a given supernova are 
meaningful but direct comparison of x 1 best-fit values for dif- 
ferent supernovae should be done only in the case of compara- 
ble signal-to-noise. 

4.3. Host galaxy subtraction 

Host galaxy subtraction is clearly an important step in the su- 
pernova spectrum identification. In our sample, some of the su- 
pernova spectra are indeed deeply contaminated by the emis- 
sion of their host. This is the case for SN 1999dt, SN 1999du 
(WHT spectrum), SN 1999dw, SN 1999dx, SN 1999dz or 
SN 1999gx. Identification of the host type is done on the full 
(i.e. SN+host) extracted spectrum by direct examination of sev- 
eral spectral features allowing to discriminate, at least broadly, 
between types. We adopt a division into three main morpho- 
logical classes, similarly to lSullivan et alJ J2003h : type corre- 
sponds to spheroids (E/SO/bulge), type 1 to early-type spirals 
(Sa/Sb) and type 2 to late-type spirals and starbursts (Sc/Stb). 



(We do not have any irregular galaxy template in our database.) 
Spectral features used for the identification (as well as for red- 
shift determination; see §4.1) include Call H&K absorption 
lines at 3934 and 3968 A, the 4000 A break (B4000), Hydrogen 
Balmer lines (mostly H/3 and Hy given the redshift range of the 
present spectra), oxygen forbidden lines [Oil] and [OIII]. We 
also consider a feature around 3850 A: the presence in early- 
type spectra of a large trough, absent in later-type spectra, due 
to CN molecular absorption blueward of 3850 A in cool stars 
(we label it T3850). This feature correlates with the strength 
of Mgb absorption at ~ 5150 A (restframe), an other useful 
metallicity indicator. Using all these features allows classifica- 
tion of all host galaxies into one of the three main categories 
defined above. 

Subtraction of the host from the supernova spectrum is per- 
formed as part of the fitting procedure by SN-fii. In two cases 
(SN 1999dr and the NOT spectrum of SN 1999du), it is possi- 
ble to use the observed host spectrum to create the model to be 
fitted. Subtraction of galaxy lines is then optimal. In all other 
cases, however, it is not possible to extract separately the host. 
We then use a host galaxy template of high signal-to-noise ratio 
appropriately redshifted (see §4.2). This offers two advantages 
over a subtraction based on using the real host spectrum: first, 
as the template used has a high S/N ratio, the resulting super- 
nova S/N is essentially not degraded by the subtraction; sec- 
ond, as the extraction procedure relies on a x 2 minimization, 
the best-fit solution gives in itself a hint of the host type (the 
'best-fit host type') that can be used as an a posteriori identifi- 
cation check by comparing to the identification based on galaxy 
features. Moreover, the ratio of the galaxy to the total super- 
nova+galaxy signal is computed by SN-fit as a byproduct of 
the fitting procedure. This ratio represents the contribution of 
the galaxy in the best-fit model in the range of overlap with 
the observed spectrum. This allows us to quantify the degree of 
contamination by the host. Although limited by the imperfect 
sampling of galactic type templates in our database, it never- 
theless gives a valuable indication on the host type and degree 
of contamination of the supernova. 



5. Spectroscopic results 

5.1. Notes on individual SN la 

As an output of the fitting procedure, a host-subtracted super- 
nova spectrum is thus obtained along with the best-fit super- 
nova template and parameters. Table |2] summarizes these re- 
sults of the fitting procedure for each supernova. Columns 3 
and 4 give the best-fit supernova template along with the per- 
centage of galaxy contamination and galaxy type. Columns 5 
and 6 give the best-fit redshift zj obtained from the fitting pro- 
cedure and the host redshift Zh respectively, z/ is obtained by 
redoing a fit of the supernova spectrum obtained after a first 
fit, cleaned from any residual galactic line. Column 7 gives the 
phase associated with the best-fit supernova+galaxy spectrum. 
As mentioned above, the phase has to be considered with a typ- 
ical ±3 day uncertainty. Finally, the minimum reduced x 2 and 
number of degrees of freedom are given in columns 8 and 9. 
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Examining, for each supernova, the solution immediately 
following the best-fit solution in the list of increasing re- 
duced^ 2 , we generally find that the phase is rather stable and 
within ±3 days of the best-fit solution. For three supernovae 
(SN 1999cj, SN 1999dr and SN 1999gx), the phase differ- 
ence is larger, up to +5 days. Note that SN 1999dr is a super- 
nova well after maximum light (its spectral features evolve less 
rapidly than around maximum). The same argument may apply 
to SN 1999cj (+9 days), although the effect is not found for the 
+9 days SN 1999dv NOT spectrum. The signal-to-noise of this 
latter spectrum is however higher than for SN 1999cj and phase 
determination should be more accurate. Finally, SN 1999gx is 
a distant SN and is dominated by the host galaxy. 

The best-fit solutions are obtained in all cases for 'normal' 
SN la templates, as opposed to peculiar SN la and other types 
(Ib/c and II). Figure ^ (upper panel) shows the best-fit red- 
shift Zf (column 5 of Tabled as a function of the redshift ob- 
tained from line identifications (the host redshift Zh, column 6 
of Table |2j. Residuals are shown on the lower panel of Fig. 
^ The two quantities agree within a few percent. Error bars 
on the host redshift are shown (they are smaller than the filled 
circle symbols and do not appear on the upper panel). When 
galactic features are weak, more difference is observed (see, 
e.g., SN 1999dy) but no clear systematic effect has been iden- 
tified (see lower panel of Fig.^. We find a r.m.s. dispersion of 
o~r.m.s ~ 0.006 for the whole sample. This value can be consid- 
ered as a typical error on the redshift determination by SN-fit. 
The average redshift of the sample is < z >= 0.348. 

We present in the lower panels of Fig.[2]to Fig.^]the output 
spectrum for the 10+2 confirmed supernovae along with the 
best-fit template overlapped on top of each spectrum. Spectra 
have been rebinned on 10A bins for visual convenience. Top 
panels show the corresponding host + SN spectra, except for 
SN 1999dr and the NOT spectrum of SN 1999du for which 
separate extraction of the host galaxy was possible (the galaxy 
alone is then shown). All spectra are presented in the observer 
frame. Residual sky lines have been removed 'by hand', but 
atmospheric absorptions and galaxy emission/absorption lines 
are left. Lines used for redshift determination are labeled. The 
symbols © indicate atmospheric absorptions. 

In Fig. 1161 and!17l we present the 12 SN classified in or- 
der of increasing spectroscopic phase. Figure ^] is for pre- 
maximum spectra, whereas Fig.El s hows past-maximum spec- 
tra. All spectra are in the restframe. Residual lines resulting 
from imperfect galaxy subtraction, as well as atmospheric ab- 
sorption lines have been removed. A polynomial least-squares 
fit using a Savitsky-Golay smoothing filter of degre e 2 with a 
window width of 60 data points (Pre ss et alJll9 92) has been 
applied to the data. This filter is designed to preserve higher 
moments within the data and is well suited for supernova spec- 
tra. For visual convenience, the spectra are shifted in flux from 
one another by an arbitrary amount. The three gray vertical 
bands show the Call, SII and Sill spectral features expected 
to be found in 'normal' SN la . In addition, three vertical solid 
lines indicate the positions of Call at 3945 A, SII at 5640 A and 
Sill at 6355 A , blueshifted by 15000 km/s (Call) and 10000 
km/s (SII, Sill). These values are typical of 'normal' SN la at 



maximum ( Be netti et all2004t) and are shown as a guide to the 
eye. 

Early spectra of 'normal' SN la are usually characterized 
by the presence of intermediate mass elements such as Si, S 
and Ca. As the spectrum evolves, they are progressively re- 
placed by features due to iron-peak ions. A strong absorp- 
tion feature of Sill visible at 6150 A in spectra around max- 
imum light is p robably the most discriminant feature against 
other SN types JWheeler & Harkneslll99rtlFmppenkolll997t 
lLeibundgutl200Cl) . However, this feature is only present in two 
of our spectra (SN 1999dr and SN 1999dv) whereas for all 
other SN it is redshifted beyond the upper limit of 7500 A even 
at moderate redshift. Other characteristic features we base our 
identification on include: 

• Presence of a W-shape visible around 5500 A (restframe) 
due to SII 5640 A 

• Presence of an absorption feature due to Sill 4130 A. 
This feature can be weak in over-l uminous supernovae of 
SN 199 IT type prior to maximum JFilippenko et alJ Il992t 
iRuiz-Lapuente et all 1 19921) . It is on the contrary abnormally 
deep and distorted in under-luminous supernovae like SN 
1991bg due to the presence of Till absorptio n in the sam e 
wavelength region i lLi et alJ E00l3 iGarnavich etaD |2004 . 
Note that for Type Ib/c supernovae, this feature presents a 
plateau- shape long ward of 4100 A u p to about 4400 A 

( Filippenkg Il997l: iBranch et alJ l2002[) similarly to its ap- 
pearance in SN 1991bg-li ke SNe JLeibundgut et ail [l993; 
iTuratto et alJl996HMazzah et alll997tlLi et alJ2001cl) . 

• Presence of a strong absorption due to Call 3945 A. This 
feature is not characteristic of SN la as it exists as well in other 
SN types and under-luminous SN la . Howe ver, it is weak in 
SN 1991T-like objects prior to maximum light dLi et alJ200ld) . 

In the following, we briefly describe the extracted spectrum 
for each supernova. 

SN 1999cj: Due to the use of the blue arm, part of the rest- 
frame UV is accessible for the two April SN la . SN 1999cj 
(Fig. |3 is at Zh = 0.362 and exhibits features characteristics 
of a SN la well after maximum. The best-fit is obtain ed with 
SN 1992A, 9 days past maximum (Kirshneret al. 1993). In par- 
ticular, suppression of the peak redward of the Sill absorption 
around 4130 A (restframe) is clearly seen. 

SN 1999ck: An early SN la at z h = 0.432 with a UV part 
from ~ 3000 A restframe (Fig. [3}- Bumps and troughs are 
clearly seen in this part of the spe ctrum. Best fit obtained with 
SN 1994D -9 days iPatat et alJl99f3) . 

SN 1999dr: This is a SN la at z h = 0.178. The host spec- 
trum has been extracted separately from the SN and used as the 
galaxy component of the fitting model. The supernova spec- 
trum is consistent with a SN la (best fit with SN 1994D) at 
24 days after maximum light. As spectra have been obtained 
on September 12th, that places the date of maximum around 
August 19, 1999 so that the supernova was starting rising when 
the reference image was taken on August 9. The discovery im- 
age was obtained 23 days after the reference. At that date, the 
flux difference between the reference and the discovery image 
still permits it to be detected as a light increasing supernova. 
Even though our detection strategy has been designed to select 
rising SNe, decreasing SNe can also be found. The aspect of 
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the spectrum (Fig. 0} is typical of a supernova in the photo- 
spheric phase a few weeks after maximum, as Fe lines develop 
and create large troughs in the spectrum. This is clearly seen 
around 6000 A 5100 A restframe) where the large dip is 
due to a blend of Fell absorption lines. The bump at * 6500 A 
replaces the W-shape due to SII seen in spectra around maxi- 
mum or earlier. 

SN 1999dt: This is the fourth most distant SN of our sam- 
ple (zh = 0.437). The spectrum is strongly dominated by the 
host galaxy and the host-subtracted spectrum presented has 
thus a rather poor S/N (see Fig.|5|l. The best-fit is for an early 
SN template (SN 1994D -9 days). The three best-fit solutions 
(in terms of x 2 ) all g ive a 'normal' SN la but the third one 
differs in phase by 7-8 days with respect to the two others: the 
best-fit solution is for -9 days while the third one yields a -1 day 
spectrum template. This 'phase instability' is not surprising as 
the spectrum is strongly dominated by the host. Unfortunately, 
as most of the candidates are detected within 1.5 arcsec of the 
host center comparable to the seeing, it is not possible to extract 
separately both components. 

SN 1999du: SN 1999du has been observed at WHT 
on September 12 and at NOT on October 2. Fig. [6] and 
show the corresponding spectra. For the NOT spectrum (Fig. 
0, the wavelength upper limit adopted for the fit is 6000 
A. Indeed, the signal redward of this value is polluted by 
the presence of the second order spectrum and is practi- 
cally useless for our fitting purposes. The same remark ap- 
plies for SN 1999dv. The October 2 spectrum is best fitted 
by a SN la (SN 19 92A) 6 days after maximum whereas a 
SN la (SN 1999ee faamuv et al J 120021) 1 at -9 days matches 
the September 12 spectrum. Taking the September 12 de- 
termination at face value, this places the date of maximum 
light on September 21, while the NOT determination leads 
to September 26 for maximum. This is marginally consistent 
given the adopted +3 days of uncertainty in the spectrum phase 
determination. Correcting for the time dilation effect expe cted 
in an expanding univer se JGoldhabe r et alJfeoOlt iRiess et all 
Il997t iFolev et al1l2005l) gives the date of maximum light on 
September 23 and September 24 for the WHT and NOT spec- 
tra, respectively, a fully consistent result. The host redshift de- 
termination agrees for the two spectra, leading to Zh = 0.260. 
Note that for the NOT spectrum, the host galaxy has been ex- 
tracted separately and used to construct the best-fit model, as 
for SN 1999dr. With this procedure, the fit is performed against 
the redshift of the supernova alone, not against the redshift of 
the supernova+galaxy template model, as for the other WHT 
spectra. This leads to a difference in z/ and Zh of 0.01 (see Table 
13, the largest difference observed for our sample (along with 
SN 1999dy). This value can be considered as an upper limit on 
the precision on the best-fit redshift zj- If we do not use the 
observed galaxy and rather use a template galaxy for the fit- 
ting model, this redshift discrepancy disappears (we then find 
Zf = 0.263). The best-fit for the WHT spectrum is obtained for 
a reddened starb urst galaxy with E(B -V)=0.4, labeled as 'Stb4' 
in our database dKinnev et alJll996l) . The same solution is ob- 
tained for the NOT spectrum when a galaxy template is used 
for the fit rather than the observed spectrum. This solution in- 
dicates a dusty star-forming galaxy. Indeed, strong emission 



lines are seen in the full spectrum. Note that excluding these 
lines from the fit leads to a solution with the same SN template 
(SN 1992 A +6 days) with a slightly lower x 2 ■ 

SN 1999dv: As for SN 1999du, both WHT and NOT spec- 
tra have been obtained for SN 1999dv and are shown on Fig. [8] 
and |g] The WHT spect rum is well fitted by a SN la - SN 2003du 
jAnupama et all20 05) at -7 days, and is only weakly contami- 
nated by its host. The host redshift is Zh = 0.186, fully consis- 
tent with the determination from the NOT host spectrum. This 
latter spectrum is best fitted by a late SN la (SN 1993 A) tem- 
plate at +9 days. Using the WHT phase, the maximum date 
falls on September 20, whereas the NOT phase determination 
leads to a maximum on September 23. These dates of maxi- 
mum are consistent within the +3 day uncertainty in the phase 
determination. Note that the phase difference found between 
the two observations (16 days) is similar within one day with 
the difference obtained for SN 1999du (15 days) but is short, 
by 3 to 4 days, of the real number of days elapsed between 
the WHT and NOT observations. Correcting for time dilation 
leads to ~ 19 days difference in the observer frame. As for 
SN 1999du, this is fully consistent with the actual observing 
dates. 

SN 1999dw: The third most distant SN la of our sample at 
Zh = 0.460. The spectrum is shown in Fig. EH It is well fitted by 
an early template (SN 1999ee -4 days) and exhibits features of 
a SN la at this date. The Call absorption feature corresponding 
to A w 3945 A restframe is clearly visible. Not however that the 
spectrum seems slightly bluer than other SN la at a comparable 
phase. 

SN 1999dx: This is a heavily buried SN la at z h = 0.269 
with a host contribution of more 75 % (Fig. II 1> . The best-fit is 
obtained for a +5 days SN la template (SN 1992A). The fit is 
poor in the range 5800 - 6200 A as the template falls below 
the observed spectrum. Interestingly, this high signal-to-noise 
supernova yields the poorest best-fit reduced x 2 of our sam- 
ple (> 1.5). Indeed, the higher the quality of the spectrum (the 
lower the noise), the more the x 2 value reflects discrepancies 
with the underlying model. This illustrates how our fitting pro- 
cedure depends on the quality of the underlying model. As ex- 
plained in §4.2, this is directly connected to the phase sampling 
and the number of the spectral SN templates in our database. 

SN 1999dy: A typical SN la at maximum with a high 
signal-to-noise ratio. The redshift is zi, = 0.215 (from galaxy 
lines, z/ being significantly lower at around zt = 0.202 ) . The 
best fit is obtained with SN 1996X dSalvo et alJl200ll) and is 
very good except around A ~ 5300 A where an emission fea- 
ture is present. Around A ~ 6300 A, the spectrum falls below 
the fit, maybe due to an uncorrected O2 atmospheric absorption 
feature around this wavelength (Fig.ll2>. 

SN 1999dz: This second most distant SN la (z h = 0.486) is 
presented in Fig.[0] The fitting procedure yields a large frac- 
tion of host galaxy and an early SN la (SN 1999ee -4 days). 
The phase is rather stable when one considers solutions of im- 
mediately higher^ 2 . Due to its fairly high redshift, features 
blueward to A ~ 4000 A (restframe) are visible. Note that the 
two peaks in the fit spectrum around 6400 - 6900 A appear 
blueshifted with respect to the peaks in the observed spectrum. 
It is difficult to determine whether this discrepancy is due to 
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an intrinsic feature of the supernova. As the galaxy contribu- 
tion to the fit model is almost 80%, it might be due to galaxy 
subtraction residuals. 

SN 1999ea: This is a SN la very similar to SN 1999dz at 
a lower redshift (zh = 0.397). The best fit is obtained for a 
SN la (SN 1994D) at -8 days, again a rather stable phase (Fig. 

m. 

SN 1999gx: The most distant SN la {z h = 0.493) of 
our sample (Fig. I15> . The SN signal is rather weak and was 
not identified during the first assessment at the telescope. 
Re-reduction of the data allowed to extract this distant SN 
( Balla nd & Mouchetll2005l) whose spectrum is best fitted by a 
normal SN la template (SN 1994D +6 days). 

5.2. Properties of the SN la 

As far as x 2 values are concerned, all the spectra presented in 
this paper are best fitted with templates of normal SN la as 
opposed to peculiar SN la . This is true even for SN 1999gx, 
even though^- 2 values are very close (see Table|4j. Given their 
signal-to-noise and redshifts, it is possible to compute a few 
spectral quantities of interest to further characterize these su- 
pernovae. This is an important study to perform as little is 
known about the spectral properties of SN la at intermediate 
redshift and the possible differences with their low-z counter- 
parts. For instance, finding spectral peculiarities in samples at 
intermediate or high redshift would shed light on the e volu- 
tion of the 'peculiarity rate' with redshift (ILi et al 1200 id) . This 
might also give some hint of how samples used for the deter- 
mination of cosmological parameters might be contaminated 
by peculiar supernovae. 

5.2.1. Velocity measurements 

Even SN la classified as 'Branch normal' do not form an 
homogeneous spectral class of objects. For instance, differ- 
ences in evolution of Call, S II and Sill abs orptio n lines veloc- 
ities are found among them dBranch et al JM 988l lHatano et alJ 
200d IBenetti et alJl2004 . SN la such as 2002bo or 1984A 
for instance exhibit rather high velocities vsui or \'caii com- 
pared to other sup ernovae considered as 'Branch normal' 
( Benetti et al. 2004). This might indicate some degree of pe- 
culiarity. However, as more spectra of nearby supernovae are 
collected and studied, there is evidence that the standard clas- 
sification in 'normal', under- or over-luminous objects is far 
too simple to accu rately reflect diversity among supernovae. 
Bene tti et al 1 d2005l) suggest a classification in 'faint' SNe with 
high-Sill velocity gradients, 'normal' SNe, also with high- 
velocity gradients but with brighter absolute magnitude, and 
SNe with low velocity gradients. These three classes are re- 
lated to the progenitor explosion mechanism. In this classifica- 
tion, 'Branch-normal' SN la are found in both the 'faint' and 
'normal' classes. 

Velocity gradient measurements are unfortunately out of 
reach for our sample, even for the two SN la for which we have 
spectra on two different dates. Call, SII and Sill absorption line 
velocities can be however measured in some cases. Results are 



shown in Table [3] Velocities have been obtained by fitting the 
corresponding line feature by a gaussian profile. Typically, the 
error is 500 km/s for Ca and 200 km/s for both SII and Sill 
velocities. We compare our velocity results to the ones of local 
supernovae la by pr oducing 'velocity vs day from B maximum' 
plots following, e.g. lBenetti et alJ 120041) . FigurefTSlshows such 
a plot for th e Call H& K featur e. Data points are taken from 
figure 1 1 of IBenetti et alJ J2004I) (see references therein). The 
date from B-maximum is estimated from the best-fit phase of 
Table |2 except for 5 SN la for which g' and r' photometry are 
available (see discussion below). 

Our different measurements are consistent with values for 
'Branch normal' SN la . The Sill velocity of SN 1999dr is 
slightly high, in between the one for SN 1994D and SN 1984A. 
SN 1999dw, SN 1999dz and SN 1999gx Call H& K veloci- 
ties are close to the ones of SN 1994D, while SN 1999dx and 
SN 1999du vcaii are closer to the one of SN 2002bo. SN 1999ea 
and SN 1999ck Call velocities are low, however consistent with 
SN 1998bu velocity. Despite this diversity, there is no hint from 
velocity measurements that any of our supernovae are peculiar. 

5.2.2. Searching for peculiar SN la 

We have used spectral information to search for peculiarities in 
the SNe of our sample (search for peculiar features in the spec- 
tra, measurements of absorption line velocities) and found that 
all our SN la are consistent with being 'normal' SN la . In order 
to assess the validity of this conclusion, we also fit our spectra 
with SN-fii using only the peculiar SN la spectral templates 
available in our database. We dispose of templates for eight pe- 
culiar SNe, both under and over-luminous with a similar phase 
sampling as for 'normal' SNe la. We list the results in Table 
0] In columns 3 and 4 the best-fit template and phase obtained 
when fitting with peculiar SN la templates are given. The corre- 
sponding reduced^ 2 and number of degrees of freedom (d.o.f.) 
are presented in columns 5 and 6. For comparison, we also re- 
port the same parameters for 'normal' SN la templates from 
Table |2 (columns 7 to 9). When fitting with peculiar templates 
alone, best-fits are systematically obtained for over-luminous 
SN l a , as seen in c olumn 3 (templates for SN 199 IT are from 
iMazzali et al] dl995h: SN 1999a a from lGaravini et alJ J2004 : 
SN 2000cx from lLietal.N2001al) \ From the comparison of the 
X 2 values in columns 5 and 8, it appears that 'normal' templates 
are favored over peculiar templates in all cases. We have per- 
formed a F-test to quantify the significance of the correspond- 
ing^- 2 difference. Column 10 gives the probability that the two 
fits are equivalent. For one SN la (SN 1999gx), the 'null hy- 
pothesis' that a solution with a peculiar SN la is equivalent to a 
solution with a 'normal' one is verified at the 2cr level. A clear 
identification is thus difficult for this object. As seen in Fig. 
[Tol SN 1999gx exhibits shallower absorptions than other SN at 
comparable phases. This might hint toward a somewhat pecu- 
liar supernova, although note that its Call velocity is normal. 

Top panel of Fig. ^] shows the best-fits obtained for 
SN 1999gx for a 'normal' and a peculiar solution. For clarity, 
the 'peculiar' solution (lower curve) has been shifted in flux by 
an arbitrary amount. The same graph for SN 1999cj, the su- 
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pernovae yielding the second highest F-test probability (57%) 
after SN 1999gx, is shown on the bottom panel. In both cases, 
visual inspection of Fig.[T9lslightlv favors the 'normal' fit over 
the peculiar, although differences tend to be small. 

For other SN la , the 'null hypothesis' is significantly re- 
jected. For SN 1999dv and SN 1999dy, two supernovae for 
which host contamination is weak, the x 2 increases dramat- 
ically, which translates in very low F-test probabilities. This 
strongly supports the hypothesis that they are 'normal' SN la . 

5.2.3. SN la Phases 

We have argued that the fitting procedure used in this paper 
allows us to constrain the phase within typically +3 days. We 
can correlate this 'spectroscopic' phase with the photometric 
phase obtained from B-band photometry for five out of the 
twelve SNe of our sample which have been photometrically 
followed-up (SN 1999dr, SN 1999du, SN 1999dv, SN 1999dx 
and SN 1999dy). The photometric phases are derived by con- 
structing and fitting the light-curves for these SNe, as described 
in our forthcoming paper. Preliminary results are given in Table 
13 A ± 1 day error bar on the photometric phase is assigned to 
each supernova. This value is likely to overestimate the true er- 
ror as light-curve fitting constrains the date of maximum within 
a fraction of a day. It however accounts for possible systemat- 
ics in the photometric analysis. We find good agreement be- 
tween the photometric and spectroscopic dates of maximum 
for the subset. When a correction for time dilation is applied, 
the agreement is even better, except for SN 1999dr. As already 
mentioned, the spectral phase determination of this SN la well 
after maximum light is likely to be less precise than for other 
supernovae. Figure [20| shows the photometric phase (column 
2 of Table as a function of the spectroscopic one (column 
3). The rms dispersion is cr rms =2.0 days. This is consistent 
with the ±3 days range we adopted for the spectroscopic phase. 
This value reduces to cr rms = 1.7 day if times dilation cor- 
rected values are considered (column 5, excluding the uncertain 
SN 1999dr). 

We also compare the phases obtained by fitting pecu- 
liar templates (Table |4} with the photometric phase for the 5 
followed-up September supernovae. The dispersion increases 
from 2.0 days ('normal') to m 6.5 days (peculiar). Again, this 
tends to show that these supernovae are 'normal'. 

5.3. Host results 

Figure|2]shows WFS reference images of host galaxies of the 
twelve supernovae presented in this paper. Each vignette is a 
0.25x0.25 square-arcmin g' image, except for SN 1999cj (B- 
band image) and SN 1999ck (r' image). Pixel size is 0.33 arsec. 
The cross indicates the position where the SN has exploded. 

Table [6] summarizes the host identification for each can- 
didate (cf. Sect. 4.3). The features used for the identification 
are presented in columns 2 and 3. Specific comments for each 
galaxy are given in column 4 while in column 5, our best de- 
termination for the host type is given. In column 6 we indicate 
the confidence we have in the identification in the form of an 



index n c whose values can be (type unsure), 0.5 (average con- 
fidence), or 1 (high confidence). Low confidence in the identifi- 
cation results from a poor signal-to-noise ratio, a fact that usu- 
ally correlates with a 'high' redshift and a poor sky subtraction, 
or with a too low contamination of the supernova by the host. 
Finally, column 7 gives the 'best-fit type' given by .SAf-fit (see 
§4.3). It is clear that in most cases, this 'best-fit type' agrees 
with the direct identification. 

For the sake of completeness, we have also computed u'-g', 
B-V and g'-r' colors for the twelve host galaxies, u', g' and r' 
magnitudes have been taken from the Sloan Digital Sky Survey 
(SDSS) on-line catalog. B- V colors have been computed us- 
ing color equations given by Fukugita et al 1 dl996l) .ErB-V)has 
been computed for each galaxy from the r' Galactic extinction 
given in the SDSS catalog and an extinction correction has been 
applied to all colors. Results are shown in Table[7] 

These derived color va lues can then be co mpared to the 
galaxy colors computed bv iFrei & Gunnl \ 19941) for their four 
Hubble galaxy types (E, Sbc, Scd and Im). To connect this clas- 
sification to the one we have adopted, we associate their type E 
to our type 0, Sbc to our type 1 and Scd and Im to our type 2. We 
assume that a Sa galaxy (type 1 in our classificat ion) has colors 
falling in between their E and Sbc types. IFrei & Gunnl d 19941) 
computed galaxy colors using th e galaxy energy dist ributions 
compiled for each Hubble type bv lColeman et alJlll980l) . Color 
values are given at four different redshifts (z = 0,0.1,0.4 and 
0.6). For the purpose of comparison, we have extrapolated Frei 
and Gunn's colors at the redshifts of our hosts, for each Hubble 
type. 

Based on this color comparison, we confirm the identifica- 
tion made from spectra for eight hosts (SN 1999cj, SN 1999ck, 
SN 1999dr, SN 1999du, SN 1999dv, SN 1999dx, SN 1999dy 
and SN 1999dz) out of twelve. Two hosts (SN 1999ea and 
SN 1999gx) have colors marginally consistent with the type de- 
rived from their spectra (at least one color does not agree with 
the spectral type). The colors of SN 1999dt and SN 1999dw 
hosts are consistent with type (two colors out of three hint 
toward a 0-type host) while they have been identified as type 1 
from their spectra. 

6. Discussion 

Based on the classification results given in Table [6] two out 
of twelve hosts are identified as type 0, seven as type 1 and 
three as type 2, corresponding to a fraction of 17%, 58% and 
25% of the sample respectively. Although we are limited by 
the small statistics of our sam ple, we can co mpare our re- 
sults with the host statistics of Sullivan et al 1 J2003I) on two 
samples of low and high redshif t. The low-z sam ple is based 
on the Calan-Tolol o sample faamuv et alJll996l) and on the 
Ries s et al.l lfl999b) local sample. All supernovae in this sam- 
ple have z < 0.01. The high-redshift sample is based on the 42 
SNe of the SCP JPerlmutter et alJll999h with redshift ranging 
from z = 0.18 to z — 0.83. Our fractions are roughly consis- 
tent with the statistics of hosts in the low-z sample of Sullivan 
et al. , who find 12%(0), 56%(1) and 32%(2). In their high-z 
sample, the fraction of E/SOs levels up with early-type spirals 
whereas the late-type spiral fraction increases: 28%(0), 26%(1) 
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and 46%(2). This trend is not seen in our sample. To check 
if consistency can be found on a subset of the high-z sample 
with similar average redshift as ours, we have selected the nine 
hosts with redshift ranging from z = 0.172 and z = 0.43 with 
an average redshift of < z >~ 0.34. However, the same trend 
as for the total high-z sample is clearly seen, with even more 
type 0. As the number of SN la of our sample is small, pos- 
sible misidentification of one host galaxy could significantly 
change our numbers. If, say, we have classified as type 1 an 
actually type galaxy (SN 1999dr for example, which has a 
red u'-g' color value consistent with the one for type 0), we 
would get 25%(0), 50%(1) and 25%(2). However, this does not 
influence the fraction of late-type spirals and we would need to 
have misclassified one type 2 as, for instance, a type 1 to al- 
leviate the discrepancy. This is rather unlikely as strong emis- 
sion lines seen in type 2 spectra should be a clear signature. In 
our sample, these strong emissions only appear for SN 1999du, 
SN 1999dv and SN 1999dy. Clearly, complementary informa- 
tion as galaxy colors and a larger set of intermediate redshift 
SN la would help to conclude. 

Taking into account various biases intrinsic to super- 
nova se arches that may affect the peculiarity rate of SN la, 
iLi et al.1 feOOld) find that 36+9% of nearby SN la are peculiar. 
Contrasting with this re sult, no peculiar SN la have been fo und 
in the h igh-z samples o flRiess et alJ Jl998 ) and P erlmutter et alJ 
ill 9991) . This might indicate evolution with redshift in the 
SN la populations, even though identification of peculiar ob- 
jects at high redshift is rendered difficult due to insufficient 
signal-to-noise and potentially severe age-bias: spectral dif- 
ferences between 'normal' and over-luminous SN la tend to 
vanish after maximum light, which introduc es an 'age-bias' in 
the detection of such over-luminous SN la dLi et al.ll2001blch . 
Assuming the fraction expected in nearby surveys holds at in- 
termediate redshift, we would expect between 3 and 5 pecu- 
liar supernovae in our sample. Our analysis yields ± 1 pecu- 
liar SN la similarly to the high-z sample. However, magnitude- 
limited surveys such as ours underestimate the rate of under- 
luminous SN la (SN 1991bg-like objects) as opposed to lo- 
cal surveys which are rather distance-limited. Over-luminous 
(SN 1 99 lT-likeO objects are usually seen in dusty star-forming 
regions dLi et alj EoOlb). Although ~0.4 mag brighter than 
'normal' SN la , they suffer from heavy extinction and might 
be under-represented in our sample. Moreover, their detection 
might be affected by the age-bias. If we assume one of our 
SN la has been incorrectly classified as a 'normal' (for exam- 
ple, SN 1999cj, SN 1999dr whose spectra have been taken well 
after maximum, or possibly SN 1999gx), the peculiarity rate in 
our sample would be ~ 8 - 9 %, a value consistent with the 12 
% expected for a magnitude-limited survey with a 20 day base- 
line, with an extra R-band extinction for SN 1991T-like objects 
of 0.8 mag and an age-bias cutoff of 7 days jLiet al.l200ld) . 

Following the first results on cosmological parameters 
based on the observation of distant supernovae, systematic 
searches for SN la have been undertaken on various telescopes 
and several samples exist today. Although aiming at observ- 
ing higher redshift SN la than ours, they often contain a few 
SN la in the intermediate redshift rang e and fairly good signal- 
to-noise spectra have been obtained. iRiess et al.l (119981) early 



sample includes five SN la in the range 0.15 < z < 0.3, while 
IPerlmutter et all d 19991) ha ve two SN la around 0.17. More re- 
centlv iBarris et al.l d20Q4t) present spectroscopy for 23 super- 
novae discovered during the IfA Deep Survey, among which 
none are in the range 0.15 < z < 0.3. Lidm an et alJ (12005) 
find two SN la in this redshift range. The large scale pro- 
gram ESSENCE (Math eson et aljE005l) presents nine spectra 
of SN la between z =0.15 and z =0.3 after two years of sys- 
tematic searches. Even more recently, discovery of supernovae 
in the SDSS survey has been reported, among which s even have 

been identified as SN la in the same redshift range dSako et alJ 

■ 1 

2005). Including our spectra, the number of SN la spectroscop- 

ically observed between 0. 15 and 0.3 amounts to thirty, all with 
similar signal-to-noise. 

Identification is always based on visual detection of char- 
acteristic spectral features such as Sill 6355A and/or com- 
parison to_w£lloteeryed local SN la . This comparison is 
direct ( Riess et alJIl 998h or involves a template fitting proce - 
du re similar as ours dBarris et alJl2004t iMatheson et alJl2005h . 
In lRiess et al.l fl998). differences between templates and real 
spectrum are small and are comparable to what we obtain with 
our se t. Th e fitting result i s of po orer quality in IBarris etafl 
d2004l) and Mathes on et alJ d2005l) . however good enough for 
the sole purpose of identification of SN la features. No galaxy 
subtraction has been attempted for the ESSENCE spectra, 
which degrades the visual aspect of the result when the full 
spectrum is fitted by local SN la templates. Note that one 
out of the nine spectr a of th e intermediate-redshift subsam- 
ple of Mathe son et al.l d2005l) shows stron g similarities with 
SN 199 IT, whereas all the five SN la of the lRiess et all ill 998) 
sample are best compared to 'normal' SN la templates. These 
results are consistent with our finding of none, or possibly one, 
peculiar SN la in our intermediate redshift sample. 

7. Conclusion 

We have presented twelve spectra of SN la supernovae taken 
at the William Herschel Telescope and the Nordic Optical 
Telescope in April, September and October 1999 during a 
search for SN la at intermediate redshifts. Five of these SN have 
redshifts between 0.1 5 < z < 0.3. This set provides high signal- 
to-noise spectra in a still largely unexplored redshift range. 

The identification was based on a ^-minimization using 
a database of galaxy and SN la templates, which allow us to 
model the observed spectra in a consistent way. Determination 
of the spectroscopic phase is reliable. Comparison with 
the photometric phase derived from the light-curves of five 
SN la that have been photometrically followed-up gives an es- 
timated +2 day uncertainty when measuring the phase from the 
matching to the spectroscopic database. 

Spectral analysis shows that most of the objects found dur- 
ing this campaign are clearly spectroscopically 'normal'. The 
observed spectra are best fitted with templates of local 'nor- 
mal' SN la, even if the possibility remains that SN 1999gx is a 
peculiar over-luminous supernova. In that case, it would be the 
most distant peculiar supernova observed so far. Velocity mea- 
surements of characteristic absorption lines such as Call, SII 
and Sill are consistent with the same measurements on 'nor- 
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mal' SN la , with SN 1999dr having a slightly higher value 
than average. 

The peculiarity rate in our sample is inconsistent with the 
rates predicted in local surveys but simila r to what is observed 
in the high-z sample of lRiess et alJ (fl998 ) and P erlmutter et alJ 
But, although unlikely, possible misidentification of a 
SN 1991T-like object as a 'normal' SN la in our sample would 
reconcile observations with the predicted local rate, provided 
that those objects suffer extra extinction. 

We conclude that the physical properties of the intermedi- 
ate redshift SN la presented in this paper are very similar, as 
far as the spectroscopic analysis is concerned, to the properties 
of their low and high redshift counterparts. 

Acknowledgements. The observations described in this paper were 
primarily obtained as visiting/guest astronomers at the INT and 
WHT, operated by the Royal Greenwich Observatory at the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de 
Astrofsica de Canarias, and the Nordic Optical 2.5 m telescope. We 
thank the dedicated staffs of these observatories for their assistance in 
pursuit of this project. We also acknowledge G. Altavilla for useful 
comments on the manuscript and G. Garavini for his help in produc- 
ing Fig. 18. We thank S. Nobili for providing us with her set of spec- 
tral templates. Some of the spectral templates used in ^SA^-fit database 
have been kindly provided to us by T. Matheson and A. Filippenko. 

References 

Aldering, G., Knop, R., & Nugent, P. 2000, AJ, 1 19, 2110 
Anupama, G. C, Sahu, D. K., & Jose, J. 2005, A&A, 429, 667 
Astier, P. & Goobar, A. 1999, IAU Circ, 7258, 1 
Balland, C. & Mouchet, M. 2005, IAU Circ, 8516, 1 
Barris, B. J., Tonry, J. L., Blondin, S., et al. 2004, ApJ, 602, 
571 

Benetti, S., Cappellaro, E., Mazzali, P. A., et al. 2005, astro-ph, 
0411059 

Benetti, S., Meikle, P., Stehle, M., et al. 2004, MNRAS, 348, 
261 

Branch, D. 1987, ApJ, 316, L81 

Branch, D., Benetti, S., Kasen, D., et al. 2002, ApJ, 566, 1005 
Branch, D., Drucker, W., & Jeffery, D. J. 1988, ApJ, 330, LI 17 
Branch, D., Lacy, C. H., McCall, M. L., et al. 1983, ApJ, 270, 
123 

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 
429, 582 

Coleman, G. D., Wu, C.-C, & Weedman, D. W. 1980, ApJS, 
43, 393 

Filippenko, A. V. 1997, ARA&A, 35, 309 
Filippenko, A. V., Barth, A. J., Matheson, T., et al. 1995, ApJ, 
450, LI 1 

Filippenko, A. V., Porter, A. C, & Sargent, W. L. W. 1990, AJ, 
100, 1575 

Filippenko, A. V., Richmond, M. W., Matheson, T., et al. 1992, 
ApJ, 384, L15 

Foley, R. J., Filippenko, A. V., E., Leonard, D., et al. 2005, 

astro-ph, 0504481 
Frei, Z. & Gunn, J. E. 1994, AJ, 108, 1476 
Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 

1748 



Garavini, G., Folatelli, G, Goobar, A., et al. 2004, AJ, 128, 387 
Garnavich, P. M., Bonanos, A. Z., Rrisciunas, K., et al. 2004, 
ApJ, 613, 1120 

Goldhaber, G, Groom, D. E., Kim, A., et al. 2001, ApJ, 558, 
359 

Hamuy, M., Maza, J., Pinto, P. A., et al. 2002, AJ, 124, 417 
Hamuy, M., Phillips, M. M., Suntzeff, N. B., et al. 1996, AJ, 
112,2408 

Hamuy, M., Pinto, P. A., Maza, J., et al. 2001, ApJ, 558, 615 
Hardin, D. 1999, IAU Circ, 7207, 3 

Hardin, D., Walton, N. A., Ellis, R. S., et al. 1999, IAU Circ, 
7182, 1 

Harkness, R. P., Wheeler, J. C, Margon, B., et al. 1987, ApJ, 
317,355 

Hatano, K., Branch, D., Lentz, E. J., et al. 2000, ApJ, 543, L49 

Home, K. 1986, PASP, 98, 609 

Howarth, I. D. 1983, MNRAS, 203, 301 

Howell, D. A., Wang, L., Supernova Cosmology Project 
Collaboration, & Supernova Factory Collaboration. 2002, 
Bulletin of the American Astronomical Society, 34, 1256 

Jha, S., Garnavich, P. M., Kirshner, R. P., et al. 1999, ApJS, 
125,73 

King, D. 1985, ING Tech Note No 31 

Kinney, A. L., Calzetti, D., Bohlin, R. C, et al. 1996, ApJ, 467, 
38 

Kirshner, R. P., Jeffery, D. J., Leibundgut, B., et al. 1993, ApJ, 
415,589 

Leibundgut, B. 2000, A&A Rev., 10, 179 
Leibundgut, B., Kirshner, R. P., Filippenko, A. V., et al. 1991, 
ApJ, 371, L23 

Leibundgut, B., Kirshner, R. P., Phillips, M. M., et al. 1993, AJ, 
105, 301 

Li, W., Filippenko, A. V., Gates, E., et al. 2001a, PASP, 113, 
1178 

Li, W., Filippenko, A. V., & Riess, A. G. 2001b, ApJ, 546, 719 
Li, W., Filippenko, A. V, Treffers, R. R., et al. 2001c, ApJ, 546, 
734 

Li, W. D., Qiu, Y. L., Qiao, Q. Y., et al. 1999, AJ, 1 17, 2709 
Lidman, C, Howell, D. A., Folatelli, G, et al. 2005, A&A, 430, 
843 

Matheson, T., Blondin, S., Foley, R. J., et al. 2005, AJ, 129, 
2352 

Matheson, T., Filippenko, A. V., Ho, L. C, Barth, A. J., & 

Leonard, D. C. 2000, AJ, 120, 1499 
Matheson, T., Filippenko, A. V., Li, W., Leonard, D. C, & 

Shields, J. C. 2001, AJ, 121, 1648 
Mazzali, P. A., Chugai, N., Turatto, M., et al. 1997, MNRAS, 

284, 151 

Mazzali, P. A., Danziger, I. J., & Turatto, M. 1995, A&A, 297, 
509 

Mazzali, P. A., Lucy, L. B., Danziger, I. J., et al. 1993, A&A, 
269, 423 

McMahon, R. G, Walton, N. A., Irwin, M. J., et al. 2001, New 

Astronomy Review, 45, 97 
Meikle, W. P. S., Cumming, R. J., Geballe, T. R., et al. 1996, 

MNRAS, 281,263 
Nobili, S., Goobar, A., Knop, R., & Nugent, P. 2003, A&A, 

404, 901 



12 



Balland et al.: Spectroscopy of twelve Type la supernovae at intermediate redshifts 



Nugent, P., Kim, A., & Perlmutter, S. 2002, PASP, 1 14, 803 
Patat, R, Benetti, S., Cappellaro, E., et al. 1996, MNRAS, 278, 
111 

Patat, R, Cappellaro, E., Danziger, J., et al. 2001, ApJ, 555, 900 
Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 
517, 565 

Perlmutter, S., Gabi, S., Goldhaber, G, et al. 1997, ApJ, 483, 
565 

Phillips, M. M., Phillips, A. C., Heathcote, S. R., et al. 1987, 
PASP, 99, 592 

Phillips, M. M., Wells, L. A., Suntzeff, N. B., et al. 1992, AJ, 
103, 1632 

Press, W. H., Flannery, B. P., Teulkolsky, S. A., & Vetterling, 
W. T. 1992, Cambridge U. Press, New York, N.Y. 

Pun, C. S. J., Kirshner, R. P., Sonneborn, G, et al. 1995, ApJS, 
99, 223 

Richardson, D., Thomas, R., Casebeer, D., Branch, D., & 
Baron, E. 2002, Bulletin of the American Astronomical 
Society, 34, 1205 

Riess, A. G, Filippenko, A. V., Challis, P., et al. 1998, AJ, 1 16, 
1009 

Riess, A. G, Filippenko, A. V., Leonard, D. C, et al. 1997, AJ, 
1 14, 722 

Riess, A. G, Filippenko, A. V, Li, W., & Schmidt, B. P. 1999a, 
AJ, 118,2668 

Riess, A. G., Kirshner, R. P., Schmidt, B. P., et al. 1999b, AJ, 
117, 707 

Riess, A. G, Strolger, L., Tonry, J., et al. 2004, ApJ, 607, 665 
Ruiz-Lapuente, P., Cappellaro, E., Turatto, M., et al. 1992, ApJ, 
387, L33 

Sainton, G. 2004a, PhD thesis, Universite Claude Bernard, 
Lyon 1 

Sainton, G. 2004b, in Exploring the Universe Contents and 
Structures of the Universe, Proc. of Moriond Conf. 

Sako, M., Romani, R., Frieman, J., et al. 2005, astro-ph, 
0504455 

Salvo, M. E., Cappellaro, E., Mazzali, P. A., et al. 2001, 

MNRAS, 321,254 
Sullivan, M., Ellis, R. S., Aldering, G., et al. 2003, MNRAS, 

340, 1057 

Turatto, M., Benetti, S., Cappellaro, E., et al. 1996, MNRAS, 
283, 1 

Turatto, M., Piemonte, A., Benetti, S., et al. 1998, AJ, 116, 
2431 

Walton, N. A., Irwin, M., McMahon, R., & Lewis, J. R. 1999, 
The Newsletter of the Isaac Newton Group of Telescopes 
(ING Newsl.), issue no. 1, p. 3-5, 1, 3 

Wheeler, J. C. & Harkness, R. P. 1990, Reports of Progress in 
Physics, 53, 1467 



Balland et al.: Spectroscopy of twelve Type la supernovae at intermediate redshifts 



13 




C 

IS) 



0.Q10E- 

0.005^ 

0.000 

-0.005 
-0.010 

m 



0.2 0.3 

Host redshift z h 



0.4 



0.5 



Fig. 1. Redshift zj from SN-fit as a function of the redshift Zh of the host galaxy (top panel), and residuals (bottom panel). Error 
bars shown are for the host redshift only. 
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Fig. 2. Top: SN 1999cj+host in the observer frame. Bottom: Galaxy-subtracted SN 1999cj, rebinned for visual convenience, with 
the best-fit template (SN 1992A +9 days) overlapped. The spectrum is shown in the observer frame and is not corrected for 
atmospheric absorptions or galactic line subtraction residuals. Note the restframe UV part accessible due to the use of ISIS blue 
arm for this supernova. 
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Fig. 3. Top: WHT spectrum of SN 1999cj+host in the observer frame. Bottom: Rebinned galaxy-subtracted spectrum of 
SN 1999ck with the best-fit template (SN 1994D -9 days) overlapped. The spectrum is shown in the observer frame and is 
not corrected for atmospheric absorptions or galactic line subtraction residuals. As for SN 1999cj, the restframe UV part is 
visible due to the use of ISIS blue arm. 
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Fig. 4. Top: WHT spectrum of SN 1999dr host galaxy alone. Bottom: Rebinned galaxy-subtracted spectrum SN 1999dr with 
the best-fit template (SN 1994D + 24 days) overlapped. The spectrum is shown in the observer frame and is not corrected for 
atmospheric absorptions or galactic line subtraction residuals. 
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Fig. 5. Top: WHT spectrum of SN 1999dt+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dt with the best-fit 
template (SN 1994D - 9 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. Note the poor [Oil], Hp and H y line subtraction and the presence of strong 
atmospheric absorptions. 
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Fig. 6. Top: WHT spectrum of SN 1999du+host. Bottom: Rebinned galaxy subtracted spectrum of SN 1999du with the best-fit 
template (SN 1999ee -9 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. 
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Fig. 7. Top: NOT spectrum of SN 1999du host alone. Due to the presence of the second order in the NOT spectra, the signal 
beyond ~ 6000 A is distorted and can not be used in the identification fitting procedure. It however shows the presence of 
galaxy lines {Hp in this case) useful for redshift determination. Bottom: Rebinned galaxy subtracted spectrum of SN 1999du 
with the best-fit template (SN 1992 A +6 days) overlapped. The spectrum is shown in the observer frame and is not corrected for 
atmospheric absorptions or galactic line subtraction residuals. The spectral range of the bottom panel does not match the one of 
the top panel. 
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Fig. 8. WHT spectrum of SN 1999dv+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dv with the best-fit tem- 
plate (SN 2003du -7 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric absorp- 
tions and galactic line subtraction residuals. 
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Fig. 9. Top: NOT spectrum of SN 1999dv+host. Due to the presence of the second order in the NOT spectra, the signal beyond 
~ 6000 A is distorted and can not be used in the identification fitting procedure. It however shows the presence of galaxy 
lines useful for redshift determination. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dv with the best-fit template 
(SN 1992 A +9 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric absorptions 
or galactic line subtraction residuals. The spectral range of the bottom panel does not match the one of the top panel. 




Fig. 10. Top: WHT spectrum of SN 1999dw+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dw with the best-fit 
template (SN 1999ee -4 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. 
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Fig. 11. Top: WHT spectrum of SN 1999dx+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dx with the best-fit 
template (SN 1992 A +5 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. Note the presence of strong atmospheric absorptions and a narrow absorption- 
like feature at 5900 A in the SN 1992 A +5 days template. 



24 Balland et al.: Spectroscopy of twelve Type la supernovae at intermediate redshifts 

SN 1999cfy + Host 

jj_jjjr — I 1 1 1 1 1 1 1 1 1 — I 1 1 1 1 — I 1 1 1 1 — i 1 1 1 1 1 1 — I 1 — r 




43M SOOO 5SW 6000 6600 7000 7500 



Fig. 12. Top: WITT spectrum of SN 1999dy+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dy with the best- 
fit template (SN 1996X day) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. 




Fig. 13. Top: WHT spectrum of SN 1999dz+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999dz with the best-fit 
template (SN 1999ee -4 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. 
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Fig. 14. Top: WHT spectrum of SN 1999ea+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999ea with the best-fit 
template (SN 1994D -8 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. 
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SN I999git + Host 




Fig. 15. Top: WHT spectrum of SN 1999gx+host. Bottom: Rebinned galaxy-subtracted spectrum of SN 1999gx with the best-fit 
template (SN 1994D -4 days) overlapped. The spectrum is shown in the observer frame and is not corrected for atmospheric 
absorptions or galactic line subtraction residuals. 
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Fig. 16. Heavily smoothed spectra of the eight pre-maximum SN la of the intermediate redshift sample ordered in a sequence of 
increasing phase (logarithmic scale). Spectra are shifted by an arbitrary amount for visual convenience, and are shown into the 
restframe. Atmospheric absorptions and galaxy line subtractions have been removed before smoothing. Grey vertical bands show 
the Call, SII and Sill features found in normal SN la . Solid vertical lines show the positions of Call at 3945 A, SII at 5640 A 
and Sill at 6355 A , blueshifted by 15000 km/s (Call) and 10000 km/s (SII, Sill). These values are typical of 'normal' SN la at 
maximum and are shown as a guide to the eye. 
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Fig. 17. Same as Fig. 16 for the six past-maximum SN la (including the two NOT spectra) of the intermediate redshift sample. 
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Fig. 18. Call velocity as a function of B-band phase. Data are taken from Figure 1 1 of Benetti et al. (2004). Calcium velocities 
for the SN la presented in this paper (plotted as filled circles) are given in Table 3. 
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Fits of SN i99Sgx with normal and peculiar SN ia templates 




4SQQ SOOO 55JO0 6000 6500 7000 7500 

A 

Fig. 19. Comparison of best fits obtained for 'normal' and peculiar SN Ia templates for two low S/N ratio supernovae for which 
identification might be questionable. Top: SN 1999gx, and bottom: SN 1999cj. For clarity, peculiar solutions have been shifted 
by an arbitrary amount in y-scale. In both cases, visual inspection slightly favors the normal solution over the peculiar one. 
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Fig. 20. Photometric phase relative to spectroscopic phase for the five photometrically followed-up SN la of our sample (top 
panel), and residuals (bottom panel). Errors shown are for the photometric phase only and assume a conservative +1 day error. 
The photometric phase is from fit of the restframe B-band light-curve and the spectroscopic phase from^- 2 minimization fitting 
of observed spectra with SN la spectral templates 
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Fig. 21. WFS reference images of host galaxies. Each vignette is a 0.25x0.25 square-arcmin cutout of the full CCD, centered on 
the position where the SN exploded (indicated by a cross). All images are in g', except for SN 1999cj (B image) and SN 1999ck 
(r' image). 
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Table 1. Log of the Spectroscopic Observations 



SN 


RA 


Dec. 


UT Date 


[R] or g' 


UT Date of 


Telescope/ 


Grating/ 


Dispersion 


Slit (arcsec)/ 




(2000) 


(2000) 


of Discovery 


magt 


Spectroscopy 


Instrument 


Grism 


(A/pix) 


Par. anglej 


1999cj 


10 14 02.46 


- 14 34.0* 


Apr 17 


[19.2] 


Apr 19 


WHT/ISIS 


R158B,R158R 


1.6/2.9 


1.0/no 


1999ck 


14 08 56.59 


- 05 50.7 


Apr 16 


[21.9] 


Apr 19 


WHT/ISIS 


R158B,R158R 


1.6/2.9 


1.0/yes 


1999dr 


23 00 17.56 


-0 05 12.5 


Sep 1 


22.1 


Sep 12 


WHT/ISIS 


R158R 


2.9 


1.5/no 


1999dt 


00 45 42.29 


+ 03 22.2 


Sep 4 


23.5 


Sep 12 


WHT/ISIS 


R158R 


2.9 


1.5/no 


1999du 


01 07 05.94 


-0 07 53.8 


Sep 8 


22.8 


Sep 12 


WHT/ISIS 


R158R 


2.9 


1.5/yes 


1999du 


" 


" 


" 


" 


Oct 2 


NOT/ALFOSC 


#4 


3.0 


2.5/yes 


1999dv 


01 08 58.96 


+ 00 24.8 


Sep 9 


21.8 


Sep 13 


WHT/ISIS 


R158R 


2.9 


1.5/yes 


1999dv 


" 


" 


" 


" 


Oct 2 


NOT/ALFOSC 


#4 


3.0 


1.0/yes 


1999dw 


01 22 52.80 


- 16 20.8 


Sep 7 


24.1 


Sep 12 


WHT/ISIS 


R158R 


2.9 


1.5/no 


1999dx 


01 33 59.45 


+ 04 15.3 


Sep 8 


22.2 


Sep 14 


WHT/ISIS 


R158R 


2.9 


1.5/no 


1999dy 


01 35 49.53 


+ 08 38.3 


Sep 8 


21.7 


Sep 14 


WHT/ISIS 


R158R 


2.9 


1.5/no 


1999dz 


01 37 03.24 


+ 01 57.9 


Sep 8 


23.4 


Sep 13 & 14 


WHT/ISIS 


R158R 


2.9 


1.5/yes 


1999ea 


01 47 26.09 


- 02 07.2 


Sep 9 


23.3 


Sep 13 


WHT/ISIS 


R158R 


2.9 


1.5/no 


1999gx 


00 34 15.47 


+ 00 04 26.2 


Sep 6 


23.3 


Sep 13 


WHT/ISIS 


R158R 


2.9 


1.5/yes 



Notes: 

1 Approximate magnitude at discovery 

* Parallactic Angle 

* Note that, due to a typographic mistake, the dec. for SN 1999cj differs by 31.5" from the value given in lHardin etalj<1999h 
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Table 2. Results from SAT-fit 



SN 


Telescope 


Best fit SN 


% (host) 1 " 


z/ 


Zh 


Spectroscopic 


Reduced 


d.o.f.* 


name 




template (reference) 








phase* 


x 2 




1999cj 


WHT 


92A (Kirshner et al. 1993) 


31 (Bulge) 


0.355 


0.362 


+9 


1.00 


859 


1999ck 


WHT 


94D (Patat et al. 1996) 


48(Sa) 


0.434 


0.432 


-9 


1.17 


502 


1999dr> 


WHT 


94D (Patat et al. 1996) 


48(own galaxy) 


0.174 


0.178 


+24 


0.95 


935 


1999dt 


WHT 


94D (Patat et al. 1996) 


80(Sb) 


0.438 


0.437 


-9 


1.24 


623 


1999du 2 


WHT 


99ee (Hamuy et al. 2002) 


65(Stb4) 


0.258 


0.260 


-9 


0.87 


853 


1999dv 


WHT 


03du (Anupama et al. 2005) 


ll(Stbl) 


0.177 


0.186 


-7 


0.99 


936 


1999dw 


WHT 


99ee (Hamuy et al. 2002) 


63(Sa) 


0.464 


0.460 


-4 


0.75 


705 


1999dx 


WHT 


92A (Kirshner et al. 1993) 


77(S0) 


0.270 


0.269 


+5 


1.66 


935 


1999dy 


WHT 


96X (Salvo et al. 2001) 


23(Sb) 


0.202 


0.215 





1.14 


934 


1999dz 


WHT 


99ee (Hamuy et al. 2002) 


79(Sb) 


0.480 


0.486 


-4 


1.17 


671 


1999ea 


WHT 


94D (Patat et al. 1996) 


42(Sa) 


0.398 


0.397 


-8 


1.12 


697 


1999gx 


WHT 


94D (Patat et al. 1996) 


64(Sa) 


0.487 


0.493 


-4 


0.93 


871 


1999du' 


NOT 


92A (Kirshner et al. 1993) 


28(own galaxy) 


0.250 


0.260 


+6 


0.81 


572 


1999dv 


NOT 


92A (Kirshner et al. 1993) 


6(Sc) 


0.181 


0.186 


+9 


1.15 


637 



Notes: 

1 %(host) = J3/(a +P) x 100 (see text) 

* Typical error is ±3 days 

* Number of degrees of freedom 

1 The observed galaxy has been used for the fit 

2 Best solution when emission lines are ignored: 99ee -9 days, 39% (Sb), z/=0.261, reduced^ 2 = 0.78, d.o.f.=936 
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Table 3. SN line velocities 



SN 


Telescope 


v c a ii 


V SII 


V SiII 


1999cj 


WHT 


-12400 






1999ck 


WHT 


-14700 






1999dr 


WHT 






-10350 


1999dt 


WHT 








1999du 


WHT 








1999dv 


WHT 




-11200 


-11300 


1999dw 


WHT 


-17400 






1999dx 


WHT 


-15650 


-9850 




1999dy 


WHT 




-11500 




1999dz 


WHT 


-16900 






1999ea 


WHT 


-15750 






1999gx 


WHT 


-16200 






1999du 


NOT 


-15200 







Notes: 

" Ca H& K ,13945 in km/s; typical error is 500 km/s 
* SII i5640 in km/s; typical error is 200 km/s 
c Sill/16355 in km/s; typical error is 200 km/s 
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Table 4. Results frormSyV-fit using peculiar SN la templates 



SN 


Telescope 




Peculiar 






Normal 




F-test 
probability (%) 


Best-fit 


Phase 


x 2 


d.o.f. 


Phase 


x 2 

A 


d.o.f. 


1999cj 


WHT 


99aa' 


+ 14 


1.025 


859 


+9 


1.005 


859 


57 


1999ck 


WHT 


99aa 


-7 


1.309 


533 


-9 


1.171 


502 


1.9 


1999dr 


WHT 


OOcx 2 


+30 


1.100 


900 


+24 


0.949 


935 


2.2 10~ 3 


1999dt 


WHT 


99aa 


-7 


1.523 


667 


-9 


1.245 


623 


6 10~ 4 


1999du 


WHT 


99aa 


-7 


0.945 


894 


-9 


0.867 


853 


1.7 


1999dv 


WHT 


99aa 


-7 


1.443 


833 


-7 


0.994 


436 


2 10- 20 


1999dw 


WHT 


99aa 


-7 


0.764 


643 


-4 


0.746 


705 


52 


1999dx 


WHT 


99aa 


-3 


1.825 


935 


+5 


1.654 


935 


0.4 


1999dy 


WHT 


99aa 


-3 


2.508 


932 





1.136 


934 


< io- 20 


1999dz 


WHT 


OOcx 


-3 


1.263 


563 


-4 


1.169 


971 


0.8 


1999ea 


WHT 


99aa 


-11 


1.321 


804 


-8 


1.117 


697 


io- 2 


1999gx 


WHT 


91T 3 


-12 


0.939 


674 


-4 


0.938 


794 


97 


1999du 


NOT 


OOcx 


+ 14 


0.883 


336 


+6 


0.812 


572 


2.7 


1999dv 


NOT 


99aa 


+ 14 


1.429 


542 


+9 


1.153 


637 


4 10~ 5 



Notes: 

1 Spectral template from Garavini et al. (2004) 

2 Spectral template from Li et al (2001c) 

3 Spectral template from Mazzali et. al. (1995) 
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Table 5. Spectral and light-curve dates of maximum for the five photometrically followed-up SN la 



SN HJD(civil date)f Uncorrected maximum date$ Corrected maximum date* 







WHT 


NOT 


WHT 


NOT 


1999dr 


2451411.45 (20/8) 


19/8 




14/8 




1999du 


2451444.83 (23/9) 


21/9 


26/9 


23/9 


24/9 


1999dv 


2451440.75 (19/9) 


20/9 


23/9 


21/9 


21/10 


1999dx 


2451428.52 (07/9) 


09/9 




08/9 




1999dy 


2451438.47 (16/9) 


14/9 




14/9 





Notes: 

t Date of B-band light-curve maximum 

£ Date estimated from the date of spectral observation corrected for the best-fit spectral phase 

* Date estimated from the date of spectral observation corrected for the best-fit spectral phase. Correction for time dilation effects has been 
applied 
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Table 6. SN host spectral classification. 



SN 


T3850 1 


Other 
indicators* 


Comments 


Host Spectral 
type* 


Confidence 
index n* 


SN-fA 


1999cj 




[Oil], s. Msb 


red, likely E/SO 





0.5 


Bulge 


1999ck 




s. [Oil], w. H/3, w. [OIII] 


poor gal. signal 


1 


0.5 


Sa 


1999dr 


n. 


s. H, K; B4000; s. Mgb 


no emission 


1 


0.5 


Own gal. 


1999dt 


n. 


s. [Oil]; s. H, K; B4000; s. Hy 


Sa/Sb 


1 


1 


Sb 


1999du 




s. [OII],H/3, [OIII] 


clear ident. Sc 


2 


1 


Stb 


1999dv 


. 


s. H/?, [OIII] 


poor gal. signal 


2 


0.5 


Stb 


1999dw 




w. [Oil] 


poor S/N 


1 





Sa 


1999dx 


y- 


s. H, K; B4000; w. [Oil]; w Hy 


clear ident. E/SO 





1 


SO 


1999dy 




[Oil]; [OIII] 


faint 


2 


0.5 


Sb 


1999dz 


n. 


w. [Oil] ; H, K 


likely Sa/Sb 


1 


0.5 


Sb 


1999ea 


n. 


s. [Oil]; [OIII]; H/3; H, K 


Sa/Sb 


1 


1 


Sa 


1999gx 




[Oil]; H, K 


Sa/Sb 


1 


0.5 


Sa 



Notes: 

1 T3850=3850A CN trough present: y.=yes, n.= no; - = not identified 

* Type 0: spheroidal (E/SO); Type 1: early-type spiral (Sa/Sb); Type 2: late-type spiral (Sc/Starburst) 

* B4000=4000A break, s.=strong, w.=weak 

* High confidence: n c =l; Average confidence: n c =0.5; Low confidence: n c =0 
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Table 7. Extinction Corrected Host Galaxy SDSS Colors 



SN 


E(B-V) 


B-V 


g'-r' 


u'-g' 


1999cj 


0.035 


1.53 


1.38 


2.12 


1999ck 


0.04 


1.33 


1.16 


0.74 


1999dr 


0.06 


1.31 


1.15 


2.12 


1999dt 


0.015 


1.53 


1.38 


3.05 


1999du 


0.03 


0.72 


0.53 


0.45 


1999dv 


0.035 


0.67 


0.47 


1.02 


1999dw 


0.04 


1.67 


1.53 


3.05 


1999dx 


0.04 


1.59 


1.44 


2.69 


1999dy 


0.03 


0.88 


0.70 


0.25 


1999dz 


0.04 


1.62 


1.47 


1.39 


1999ea 


0.035 


1.03 


0.85 


1.22 


1999gx 


0.02 


1.18 


1.01 


0.86 
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Appendix A: SN-1 it spectral database 

A total of 242 spectral templates have been collected from various sources into the SN-fit database used for the analysis of the spectra presented 
in this paper. These include 167 type la spectra (both 'normal' and peculiar 2 ), 64 type Ib/c and II spectra and 11 galaxy spectra. Mos t of the 
type la and type II templates originate from the Asiago Catalog 3 and have been retrieved from the SUSPECT 4 database l Richardson et al. 
I2OO2I) with the exception of SN 1981 B, SN 1992A, as well as some spectra of SN 1990N and SN 1991T. Most typ e Ib/c templates are from 
Matheson (private communication and M atheson et alj|200ll) ). Galaxy templates are taken from iKinnevetalJi 19961) and cover a large spectral 
range (from 1200 to 10000 A). Morphological types include bulge, E, SO, Sa, Sb, Sc, and starbursts (Stb) with various amount of reddening 
ICalzetti et alll994tlKinnev et alll996l) . 

The selected supernova spectra have: i) a large wavelength coverage in the visible - typically 4000 to 8000 A restframe, ii) a good signal 
to noise, iii) a continuous energy distribution in the wavelength range considered (templates with gaps in their energy distribution have been 
discarded). All spectra have been de-redshifted and de-reddened when necessary. The spectra phases cover the range -15,45 days - with respect 
to maximum light - as uniformly as possible. A few later phase spectra (phase greater than 45 days) have been included to the database for 
peculiar SN la (4 spectra) and Ib/Ic/II (8 spectra). Table lA. f1 and lA.2l show the supernova templates available in .SJV-fit database along with the 
corresponding number of spectra, for type la and Ib/Ic/II respectively. 

The phase coverage of SN la and SN Ib/Ic/II templates available in the database is shown as histograms in Fig. IA.ll and lA.2l The phase 
ranges from -15 to +45 days. Note however that the total number of spectra indicated in the figures takes into account the spectral templates 
with phases greater than 45 days, which are not shown here. 

In Fig. lA. II both the phase distribution for 'normal' SN la (top panel) and peculiar SN la (91bg-like SN la : middle panel; 91T-like SN la : 
bottom panel) are shown. Phases for normal SN la templates are uniformly spread between -15 and 40 days, except for a gap of three days after 
day 20. This gap will be filled in the future as more spectra of low-redshift supernovae become available. Note also the absence of spectra at 
day -12. 

91bg-like SN la templates are found mostly around maximum. The sampling is smooth between -5 and +8 days, with the exception of day 
-2 for which templates are missing. No template is available between day 9 and day 15. The distribution of 91T-like SN la templates peaks a 
few days before maximum. Note the absence of spectra at days 2, 3 and 4. After day +11, the sampling is scarcer with a phase resolution of 2 
days up to about day 30. 

For completeness, we also show in Fig. lA.2l the phase distribution for type Ib/Ic and II. It can be seen that spectral phases for these types 
are best represented in the database between -5 and +15 days. 

Note that SN-fit database is in constant evolution. New spectra are regularly included as they become available to the community. 



2 Peculiar la include 91bg-like and 91T-like SN la 

3 http://web.pd.astro.it/supern/ 

4 http : //bruf ord . nhn . ou . edu/~suspect/ 
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Table A.l. SN la template spectra in .SJV-Fit Database 



SN 


Type 


Number of 


Reference 


Name 




spectra 




SN 1981B 


la 


6 


Branch et al. (1983) 


SN 1986G 


91bg-like 


5 


Phillips etal. (1981) 


SN 1990N 


la 


7 


Leibundsut et al. (1991): Phillips et al. (1992): Mazzali et al. f 1993") 


SN 1991T 


91T-like 


9 


Filippenko et al. (1992) 


SN 1991bg 


91bg-like 


6 


Phillips et al. (1992): Mazzali et al. (1995) 


SN 1992A 


la 


12 


Kirshner et al. ^1993) 


SN 1994D 


la 


18 


Meikle et al. (1996): Patat et al. (1996) 


SN 1996X 


la 


6 


Salvo etal. (2001) 


SN 1997br 


91T-like 


10 


Li et al. (1999) 


SN 1997cn 


91bg-like 


1 


Turatto et al. (1998) 


SN 1998bu 


la 


12 


Jhaetal. (1999) 


SN 1999aa 


91T-like 


10 


Garavini et al. (2004) 


SN 1999by 


91bg-like 


14 


Garnavich et al. (2004) 


SN 1999ee 


la 


11 


Hamuv et al. (2002) 


SN 2000cx 


91T-like 


24 


Li et al. (2001a) 


SN 2002bo 


la 


12 


Benetti et al. (2004) 


SN 2003du 


la 


4 


Anunama et al. (2005) 


Total peculiar la = 79 


Total 'normal la' = 88 


Total = 167 
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Table A.2. SN lb, Ic and II template spectra in SN-Fit Database 
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SN 


Type 


Number of 


Reference 


Name 




spectra 




SN 1983N 


lb 


2 


Harkness et al. (1987) 


SN 1984L 


lb 


5 


Harkness et al. (1987) 


SN 1987A 


ILpec 


4 


Pun et al. (1995) 


SN 1987M 


Ic 


1 


FiliDDenko et al. (1990) 


SN 1990B 


Ic 


4 


Matheson et al. (2001) 


SN 1990U 


Ic 


3 


Matheson et al. (2001) 


SN 1993 J 


lib 


6 


Matheson et al. (2000) 


SN 19941 


Ic 


7 


FiliDDenko et al. (1995) 


SN 1997ef 


Ibc_pec 


3 


Matheson et al. (2001) 


SN 1998bw 


Ic 


7 


Patatetal. (2001) 


SN 1998dt 


lb 


2 


Matheson et al. (2001) 


SN 1999di 


lb 


2 


Matheson et al. (200 1 ) 


SN 1999dn 


lb 


3 


Matheson etal. (2001) 


SN 1999em 


Hp 


11 


Hamuvetal. (2001) 


SN 1999ex 


Ibc 


3 


Hamuv et al. (2002) 


SN 2000h 


lb 


1 


Branch et al. (2002) 


Total Ib/Ic = 43 


Total II =21 


Total Ib/Ic/II = 64 
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Fig.A.l. Phase distribution of type la templates in SN-Ht database. Up: Normal SN la . Middle: 91bg-like SN la . Bottom: 
91T-likeSNIa. 
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Phase histogram for all Ib/c and II 
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Fig. A.2. Phase distribution of type lb, Ic and II templates in SN-Ht database. 



